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Abstract  Spin-orbit  coupling  is  often  described  in  an  approach  known  as  “the  MacDonald 
torque ”,  which  has  long  become  the  textbook  standard  due  to  its  apparent  simplicity.  Within 
this  method,  a  concise  expression  for  the  additional  tidal  potential,  derived  by  MacDonald 
(Rev  Geophys  2:467-541,  1964),  is  combined  with  a  convenient  assumption  that  the  quality 
factor  Q  is  frequency-independent  (or,  equivalently,  that  the  geometric  lag  angle  is  constant 
in  time).  This  makes  the  treatment  unphysical  because  MacDonald’s  derivation  of  the  said 
formula  was,  very  implicitly,  based  on  keeping  the  time  lag  frequency-independent,  which 
is  equivalent  to  setting  Q  scale  as  the  inverse  tidal  frequency.  This  contradiction  requires 
the  entire  MacDonald  treatment  of  both  non-resonant  and  resonant  rotation  to  be  rewritten. 
The  non-resonant  case  was  reconsidered  by  Efroimsky  and  Williams  (Cel  Mech  Dyn  Astron 
104:257-289,  2009),  in  application  to  spin  modes  distant  from  the  major  commensurabil- 
ities.  In  the  current  paper,  we  continue  this  work  by  introducing  the  necessary  alterations 
into  the  MacDonald-torque-based  model  of  falling  into  a  1-to-l  resonance.  (The  original 
version  of  this  model  was  offered  by  Goldreich  (Astron  J  71 : 1-7,  1966).  Although  the  Mac¬ 
Donald  torque,  both  in  its  original  formulation  and  in  its  corrected  version,  is  incompatible 
with  realistic  rheologies  of  minerals  and  mantles,  it  remains  a  useful  toy  model,  which 
enables  one  to  obtain,  in  some  situations,  qualitatively  meaningful  results  without  resort¬ 
ing  to  the  more  rigorous  (and  complicated)  theory  of  Darwin  and  Kaula.  We  first  address 
this  simplified  model  in  application  to  an  oblate  primary  body,  with  tides  raised  on  it  by 
an  orbiting  zero-inclination  secondary.  (Here  the  role  of  the  tidally-perturbed  primary  can 
be  played  by  a  satellite,  the  perturbing  secondary  being  its  host  planet.  A  planet  may  as 
well  be  the  perturbed  primary,  its  host  star  acting  as  the  tide-raising  secondary).  We  then 
extend  the  model  to  a  triaxial  primary  body  experiencing  both  a  tidal  and  a  permanent- 
figure  torque  exerted  by  an  orbiting  secondary.  We  consider  the  effect  of  the  triaxiality  on 
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both  circulating  and  librating  rotation  near  the  synchronous  state.  Circulating  rotation  may 
evolve  toward  the  libration  region  or  toward  a  spin  faster  than  synchronous  (the  so-called 
pseudosynchronous  spin).  Which  behaviour  depends  on  the  orbit  eccentricity,  the  triaxial 
figure  of  the  primary,  and  the  mass  ratio  of  the  secondary  and  primary  bodies.  The  spin  evo¬ 
lution  will  always  stall  for  the  oblate  case.  For  libration  with  a  small  amplitude,  expressions 
are  derived  for  the  libration  frequency,  damping  rate,  and  average  orientation.  Importantly, 
the  stability  of  pseudosynchronous  spin  hinges  upon  the  dissipation  model.  Makarove  and 
Efroimsky  (Astrophys  J,  2012)  have  found  that  a  more  realistic  tidal  dissipation  model  than 
the  corrected  MacDonald  torque  makes  pseudosynchronous  spin  unstable.  Besides,  for  a 
sufficiently  large  triaxiality,  pseudosynchronism  is  impossible,  no  matter  what  dissipation 
model  is  used. 

Keywords  Bodily  tides  •  Land  tides  •  Moon  •  Libration  ■  Spin-orbit  resonance  • 

Planetary  satellites 

1  Motivation 

Bodily  tides  in  a  near-spherical  homogeneous  primary  perturbed  by  a  point-like  secondary 
are  described  by  the  theory  developed  mainly  by  Darwin  (1879,  1880)  and  Kaula  (1966). 
Sometimes  this  theory  is  referred  to  as  the  Darwin  torque.  Based  on  a  Fourier-like  expansions 
of  the  perturbing  potential  and  of  the  tidally-induced  potential  of  the  disturbed  primary,  their 
theory  permits  an  arbitrary  rheology  of  the  primary. 

Sometimes  a  much  simpler  empirical  model,  offered  by  MacDonald  (1964)  and  often 
named  as  the  MacDonald  torque,  is  used  in  the  literature  for  obtaining  very  approximate 
but  still  qualitatively  reasonable  description  of  tidal  evolution.  This  model  can  be  derived 
from  the  Darwin-Kaula  theory,  under  several  simplifying  assumptions.  Among  these  is  the 
assumption  that  the  tidal  quality  factor  Q  of  the  perturbed  primary  should  scale  as  the  inverse 
of  the  tidal  frequency.  Being  incompatible  with  the  rheologies  of  actual  minerals,  this  key 
assumption  prohibits  the  use  of  the  MacDonald  model  in  long-term  orbital  calculations. 

Despite  the  unrealistic  rheology  instilled  into  the  MacDonald  model,  the  model  remains  a 
lab  which  permits  one  to  gain  qualitative  understanding  of  tidal  evolution  (rotational  and 
orbital),  without  resorting  to  the  lengthy  calculations  required  in  the  accurate  Darwin- 
Kaula  approach.  It  should  be  noted  however  that  employment  of  the  MacDonald  model  needs 
some  care.  Historically,  the  orbital  calculations  performed  with  aid  of  this  model  by  one  of 
its  creators,  MacDonald  (1964),  contained  an  inherent  contradiction.  MacDonald  began  his 
paper  with  deriving  a  concise  expression  for  the  additional  tidal  potential,  and  then  com¬ 
bined  this  expression  with  a  convenient  assumption  that  the  geometric  lag  angle  is  constant 
(or,  equivalently,  that  the  tidal  quality  factor  is  a  frequency  independent  constant).  This 
made  the  treatment  unphysical  because  MacDonald’s  derivation  of  the  said  formula  was, 
very  implicitly,  based  on  keeping  the  time  lag  frequency-independent,  which  is  equivalent  to 
setting  Q  to  scale  as  the  inverse  tidal  frequency.  This  contradiction  made  his  theory  incon¬ 
sistent.  The  said  oversight  has  then  been  repeated  many  a  time  in  the  literature  (Kaula  1968, 
Eq.  4.5.37;  Murray  and  Dermott  1999,  Eq.  5. 14). 1  In  particular,  the  MacDonald  method  was 
used  by  Goldreich  (1966)  in  his  theory  of  dynamical  evolution  near  the  1:1  spin-orbit  reso- 

1  Due  to  an  error  in  our  translation  from  German,  we  mis-assumed  in  our  previous  papers  (Efroimsky  and 
Williams  2009;  Efroimsky  2012a;  that  Gerstenkorn  1955)  had  based  his  development  on  a  constant- (9  model. 
Therefore  we  stated  that  his  theory  contained  the  same  genuine  inconsistency  as  the  theory  by  MacDonald 
(1964).  Accurate  translation  of  the  work  by  Gerstenkorn  (1955)  has  shown  that  his  method  was  based  on 
a  constant-time-lag  model.  Therefore  we  retract  our  statement  about  Gerstenkorn’s  approach  sharing  the 
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nance.  We  reconsider  this  theory,  employing  the  corrected  version  of  the  MacDonald  torque, 
i.e.,  setting  the  quality  factor  to  scale  as  inverse  frequency. 

2  Linear  bodily  tide  in  a  near-spherical  primary 

Consider  a  near-spherical  primary  of  radius  R  and  a  secondary  of  mass  M*ec  located  at 
r  =  (r* ,  (j>* ,  A*),  where  r*  >  R.  The  tidal  potential  created  by  the  secondary  alters  the 
primary’s  shape  and,  as  a  result,  its  potential.  For  linear  tides,  the  amendment  to  the  primary’s 
exterior  potential  is  known  (e.g.,  Efroimsky  and  Williams  2009;  Efroimsky  2012a, b)  to  read 
as 

°o  Rll+1  1  (l— mV. 

U  (?)  =  -  GM*ec  ^  ki  I+]  l+]  X  ;(2  -  80m)Plm  (sin  c/>)Plm(sin<P*) 

1=2  r  r*  m=0  '  '  m'' 

cos/n(A  — A*),  (1) 

8jj  being  the  Kronecker  delta,  G  =  6.7  x  10~u  m3  kg^3s^2  being  Newton’s  gravity  constant, 
and  y  being  the  angle  between  the  vectors  r  *  and  r  pointing  from  the  primary’s  centre.  As 
agreed  above,  r  =  (r* ,  (f>* ,  A*)  denotes  the  position  of  the  perturber,  while  r  =  (r,tp.  A) 
is  an  exterior  point,  at  a  radius  r  >  R,  where  the  tidal  potential  amendment  U ( r )  is  mea¬ 
sured.  The  longitudes  A,  A*  are  reckoned  from  a  fixed  meridian  on  the  primary,  while  the 
latitudes  <f>,  <f>*  are  reckoned  from  its  equator.  Indices  /  and  m  are  traditionally  referred  to 
as  the  degree  and  order,  accordingly.  The  associated  Legendre  functions  Pim{x)  (termed 
associated  Legendre  polynomials  when  their  argument  is  sine  or  cosine  of  some  angle)  are 
introduced  as  in  Kaula  (1966,  1968),  and  may  be  called  unnormalised  associated  Legendre 
functions,  to  distinguish  them  from  their  normalised  counterparts.2  The  Love  numbers  lq  can 
be  calculated  from  the  geophysical  properties  of  the  primary. 

A  different  formula  for  the  tidal-response-generated  change  in  the  potential  was  suggested 
by  Kaula  (1961,  1964).  Kaula  devised  a  method  of  switching  variables,  from  the  spherical 
coordinates  to  the  orbital  elements  (a*,  e* ,  i* ,  £2*,  of ,  A4*  )  and  (  a,  e,  i,  Q,  a>,  A4)  of 
the  secondaries  located  at  r  and  r.  The  goal  was  to  explore  how  a  tide-raising  secondary 
at  r  acts  on  a  secondary  at  r  through  the  medium  of  the  tides  it  creates  on  their  mutual 
primary. 

The  development  enabled  Kaula  to  process  (1)  into  a  series,  which  was  a  disguised  form 
of  a  Lourier  expansion  of  the  tide.  Interestingly,  Kaula  himself  never  referred  to  that  expan¬ 
sion  as  a  Fourier  series,  nor  did  he  ever  write  down  explicitly  the  expressions  for  the  Fourier 
modes.  At  the  same  time,  the  way  in  which  Kaula  introduced  the  phase  lags  indicates  that 
he  was  aware  of  how  the  modes  were  expressed  via  the  perturber’s  orbital  elements  and  the 


Footnote  1  continued 

inconsistency  of  MacDonald’s  theory.  We  also  thank  Hauke  Hussmann  and  Peter  Noerdlinger  for  their  kind 
help  in  translating  excerpts  from  Gerstenkom’s  work. 


The  Legendre  polynomials  may  be  defined  through  the  Rodriguez  formula  P\{x)  =  yx  —  I  J  . 

In  most  literature,  the  associated  Legendre  functions  are  introduced,  for  a  nonnegative  m ,  as 


Pin 


/  ~\m/ 2  rf'"  /  -,\m/z  rtf'" 

(*>  =  (‘  -  *  )  ^  and  =  ("U  {l~x)  PlW' 


mil  dm 
dx " 


so  that  P/m  ( x )  =  (— l)m  Pjn  (*).  The  above  definition  agrees  with  the  one  offered  by  Abramowitz  and  Stegun 
(1972,  p.  332).  A  different  convention  is  accepted  in  those  books  (e.g.,  Arfken  and  Weber  1995,  p.  623)  where 
P^ix)  lacks  the  (— l)m  factor  and  thus  coincides  with  Pim(x). 
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primary’s  rotation  rate.  The  original  works  by  Kaula  (1961,  1964)  were  written  in  a  terse 
manner,  with  many  technicalities  omitted.  A  comprehensive  elucidation  of  his  approach  can 
be  found  in  the  Efroimsky  and  Makarov  (2012).  Referring  the  reader  to  that  paper  for  details, 
here  we  cite  only  the  resulting  formula  for  the  secular  part  of  U,  in  the  special  case  when  the 
tide-raising  secondary  itself  experiences  perturbation  from  the  tide  it  creates  on  the  primary 
(so  r  *  =  r): 


U{sec\r) 


G  MSec 
a 


l 

Z 

m=  0 


(/  —  m)\ 
(l  +  m)\ 


(2 


So  m) 


l  oo 

x2X,(o  £  kl  ((Oj  rn  ,JC/ )  COS  €[  ((O/fjipq  )  , 

p= 0  <7=-oo 


(2) 


where  /,  m,  p,  q  are  integers,  Fimp(i )  are  the  inclination  functions  (Gooding  and  Wagner 
2008),  Gipq(e)  are  the  eccentricity  polynomials  coinciding  with  the  Hansen  coefficients 

(-1-1),  (l-2p) 

X lt  2p  (  )  (e),  while  the  superscript  “sec”  means:  secular. 

The  dynamical  Love  numbers  k /  and  the  phase  lags  ej  are  functions  of  the  Fourier  modes 


coimpq  =  (/  -  2 p)  d>  +  {l  -  2p  +  q)  M  +  m  {£l  -  6),  (3) 

9  being  the  primary’s  sidereal  angle,  and  9  being  its  angular  velocity.  While  the  tidal  modes 
(3)  can  be  of  either  sign,  the  physical  forcing  frequencies 


Xlmpq  =  I  coimpq  I  =  I  (/  -  2 p)  o>  +  (Z  -  2p  +  q)  M  +  m  (f2  -  9)  \  (4) 


at  which  the  stress  and  strain  oscillate  in  the  primary  are  positive  definite. 

A  partial  sum  of  series  (2),  with  |/|,  \q\,  |j|  <  2,  was  offered  earlier  by  Darwin  (1879). 
An  explanation  of  Darwin’s  method  in  modern  notation  can  be  found  in  Ferraz-Mello  et  al. 
(2008). 3 

The  power  of  the  Darwin-Kaula  approach  lies  in  its  compatibility  with  any  rheology,  i.e., 
with  an  arbitrary  form  of  the  mode-dependence  of  the  products  k\  cos  e; .  It  can  be  dem¬ 
onstrated  that,  for  a  homogeneous  near-spherical  primary,  the  functional  form  of  the  mode 
dependence  of  the  product  k\  cos  eimpq  is  defined  by  index  l  solely:  k\  ( a>impg )  cos  e/  {u>impq), 
the  other  three  indices  being  attributed  to  the  tidal  mode.  One  can  assume  that  this  prod¬ 
uct  depends  not  on  the  tidal  mode  cDi,npq  but  on  the  positive  definite  physical  frequency 
Xlmpq  ■  This  however  will  require  some  care  in  derivation  of  the  tidal  torque  from  the  above 
expressions  for  the  potential — see  Efroimsky  (2012a,b)  for  details. 

We  had  to  write  down  the  secular  part  of  the  Kaula  expansion  of  tide,  because  we  shall 
use  it  as  a  benchmark  wherewith  to  compare  the  empirical  expression  by  MacDonald  (1964), 
which  we  shall  derive  below  in  an  accurate  manner. 

Lacking  the  ability  to  accommodate  an  arbitrary  rheology,  the  MacDonald  approach 
(or  the  MacDonald  torque)  produces,  after  a  necessary  correction,  a  simple  method  which 
can  sometimes  be  employed  for  getting  qualitative  understanding  of  the  picture. 


3  Tidal  torque 

Consider  a  secondary  body  of  mass  Msec,  located  relative  to  its  primary  at  r  =  (r,  k,  <f>), 
where  4>  is  the  latitude,  and  k  denotes  the  longitude  reckoned  from  a  meridian  fixed  on 

'  In  Ferraz-Mello  et  al.  (2008),  the  meaning  of  notations  r  and  r  *  is  opposite  to  ours. 
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the  primary.  Let  U  stand  for  the  tidal-response  amendment  to  the  primary’s  potential.  This 
amendment  can  be  generated  either  by  this  secondary  itself  or  by  some  other  secondary  of 
mass  M*ec  located  at  r*  =  (r* ,  X*,  <p*).  In  either  case,  the  primary’s  tidal  response  to  the 
gravity  of  the  secondary  will  render  a  tidal  force  and  torque  acting  on  the  secondary  of  mass 
Msec.  The  torque’s  component  perpendicular  to  the  equator  of  the  primary  will  be  given  by: 


Tz 


MSec 


dU(r) 

dX 


(5) 


We  would  reiterate  that  (5)  is  a  component  of  the  torque  wherewith  the  primary  acts  on  the 
secondary  of  mass  Msec. 

Then  its  negative  will  be  the  appropriate  (i.e.,  orthogonal  to  the  primary’s  equator)  com¬ 
ponent  of  the  torque  wherewith  the  secondary  acts  on  the  primary: 


Tz(r)  =  -  Tz 


MSec 


dU(r) 

dX 


(6) 


Derivation  of  formulae  (5-6)  is  presented  in  Appendix  A.  These  expressions  are  convenient 
when  the  tidal-response  potential  amendment  U  is  expressed  through  the  spherical  coordi¬ 
nates  r,  X,  cp  and  r*  ,X*  ,cp* ,  as  in  formula  (1). 

Whenever  the  tidal  response  is  expressed  as  a  function  of  the  orbital  elements  of  the 
secondary  and  the  sidereal  angle  9  of  the  primary,  it  is  practical  to  write  down  that  the 
perpendicular-to-equator  component  of  the  torque  acting  on  the  primary  as 


Tz{r) 


-  Ms, 


dU(r) 

de 


(7) 


9  standing  for  the  primary’s  sidereal  angle  (Efroimsky  2012a,b). 

We  prefer  to  employ  the  terms  primary  and  secondary  rather  than  planet  and  satellite. 
This  choice  of  terms  is  dictated  by  our  intention  to  apply  the  below-developed  machinery 
to  research  of  tidal  dissipation  and  spin  evolution  of  a  satellite.  In  this  setting,  the  satellite 
is  effectively  playing  the  role  of  a  tidally-distorted  primary,  its  host  planet  acting  as  a  tide¬ 
generating  secondary.  Whenever  the  method  is  applied  to  exploring  the  problem  of  planet 
despinning,  the  planet  is  understood  as  a  primary,  the  host  star  being  a  tide-raising  secondary. 


4  MacDonald  (1964) 

When  it  comes  to  taking  dissipation  into  account,  expression  (1)  turns  out  to  be  far  more 
restrictive  than  (2),  in  that  (1)  becomes  applicable  to  a  very  specific  rheological  model.  This 
happens  because  a  straightforward4  option  of  instilling  the  delay  into  (1)  is  to  replace  in 
this  expression  the  perturber’s  coordinates  r*{t) ,  <p*{t) ,  X*(t)  with  their  delayed  values, 
r*(t  —  At) ,  (p*(t  —  At)  ,  X*(t  —  At).  For  example,  instead  of  cos m(X  —  X*)  we  should 
employ 

4  A  more  accurate  treatment,  which  cannot  be  employed  within  the  MacDonald  model  but  is  implementable 
within  the  Darwin-Kaula  approach,  is  to  expand  the  tide-raising  potential  W  and  the  tidal-response  potential 
change  U  into  Fourier  modes  n>impq ,  and  then  to  introduce  the  Love  numbers  kjmptl  =  k\  {<oimpq ) ,  phase  lags 
€lmpq  —  e/  and  time  lags  Atimpq  =  At/(o);mpfy).  This  formalism  (explained  in  detail  in  Efroimsky 

2012a,b)  enables  one  to  express  the  so-introduced  Love  numbers  through  the  rheological  properties  of  the 
primary's  mantle,  and  thereby  to  model  adequately  the  frequency-dependence  of  these  Love  numbers. 

An  intermediate,  purely  empirical,  option  would  be  to  introduce  “Love  numbers”  k[m,  as  if  they  were  func¬ 
tions  of  both  the  degree  /  and  order  m .  This  idea  is  implemented  in  the  IERS  Conventions  on  the  Earth 
rotation  (Petit  and  Luzum  2010).  In  the  LLR  (Lunar  Laser  Ranging)  integration  software,  tides  in  the  Earth 
are  parameterised  by  k[m  and  A (/,„  with  1  =  2  and  m  =  0, 1,  2  (Standish  and  Williams  2012). 
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( delayed ) 


cos  I  m  X  —  m  X 


=  cos  (  m  X  —  (mX  —  m  X  A t  ) 


(8) 


insertion  whereof  into  (1)  yields: 

k,GM* 


U(r)  =  ~X 


21+1  / 

R  (/  —  mV- 

2_,  773 - 77  (2  ~S0 m)  Pin,  ( Sin  0(f)) 


~2  Ht)l+1  ti +l  +  m)[ 

P/,„(sin0*(t  —  At))  cosm(X  —  X*  +  X*  At).  (9) 

Expressing  the  longitude  via  the  true  anomaly  v, 

/.  —  —  Q-\-Q-\-co~\-v-\-  0{i2)  =  —  6  T  £2  4-  T  AA  -t-  2  e  sin  Ad 

+  0(e2)  +  0(i2),  (10) 


and  neglecting  the  apsidal  and  nodal  precessions,  we  obtain: 


cos  (  (m  X  —  mX  )  +  mX  At 

■(< 


) 


=  cos  ^  (m  X  —  mX  )  +  m  (v*  —  6*)  At  +  O  ( i *2)  ^  (11) 

or,  in  terms  of  the  mean  anomaly: 

cos  ^(/?iA  —  mX  )  +  m  X*  Af^  =cos  ^(mX  —  mX  )+m(n*  —  6*)At  +  2me*n*  At  cosA/f* 


+  0(e*2)  +  O  (i*2)+  0(i2))  . 


(12) 


This  enables  us  to  write  down  the  potential  as 


U(r)  =  —GM*ec  ]>> 


R 


2  l+l 


■r— i  (Z  —  m)\ 

/  ,  „  .  ,,  (2  -  S0m)  P/m  (sin  0(f)  ) 


Ht)1+1  r\t  -  At)l+1  ^o(/  +  m)! 

P/m  (sin</>*  (f  —  At))  cos  (X  —  A.  )  +  /n  (v*  —  6*)  At  +  0(i2)  +  0(i*2)^ 


(13) 


where  P/,„  (sin  </>(?)  )P/,„(sin0*(t  —  At) )  may  be  replaced,  in  the  order  of  0(i 2)  +  0(i*2)  + 
0  (ii* )  ,  with  Pi,n  (0)  Pi,, i  (0)  .  A  further  simplification  can  be  achieved  by  taking  into  account 
the  /  =  2  contribution  only.  Omitting  the  /.- independent  term  with  m  =  0  (in  the  expres¬ 
sion  for  the  torque,  m  will  become  a  multiplier  after  the  differentiation  of  U  (?)  with  respect 
to  X),  and  omitting  the  m  =  1  term  (as  P21  (0)  =  0  ),  we  arrive  at  the  expression 

U(r)  =  --  GI^sec  kl  R  [l  +  0(i2)  +  0(i*2)  +  0(ii*)  1  cos  ((2X  -  2X*  ) 

4  r(0  r*(t-  At)  1  J  V 

+  2  (v*  -  0*)  At  +  0(i2)  +  0(i*2)  )  .  (14) 

In  the  special  case  when  the  tide-raising  satellite  is  the  same  body  as  the  tidally-perturbed  one 
(i.e.,  when  r(f)  =  r*(t) ,  Msec  =  M*ec,  and  X  =  X*),  this  expression  happens  to  coincide, 
in  the  leading  order  of  i,  and  i*.  with  the  potential  employed  by  MacDonald  (1964).  Thus 
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we  have  reproduced  his  empirical  approach,  by  starting  with  the  rigorous  formula  (1),  and 
by  performing  the  following  sequence  of  approximations: 

•  First,  when  accommodating  dissipation,  we  set  all  the  time  delays  to  be  equal,  for  all  the 
physical  frequencies  Ximpq  —  I  wimpq  I  involved  in  the  tide: 

Afl,nPq  =  A«Xlmpq)  =  At ■  (15) 

This  point  is  explained  in  great  detail  in  Efroimsky  and  Makarov  (2012). 

•  Second,  we  assume  the  smallness  of  the  inclinations  and  lag,  through  the  neglect  of  the 
relative  errors  0(;'2),  0(i*2),  and  0(ii*). 

•  Third,  we  truncate  the  series  by  leaving  only  the  l  =  m  =  2  term. 


These  three  steps  take  us  from  (1)  to  the  MacDonald  approximation 

3  GM*  k2  R  (  *  \ 

U (?)  «  -  -  - ,  sec  - T  cos  (2k  -  2k  +  e  ) 

^  r(t)  r*(t  —  At) 


(16a) 


where  r*(t  —  At)  can  be  replaced,5  in  the  order  of  0(e*n*  At),  with  r*(t).  With  the  latter 
simplification  implemented,  and  with  the  tide-raising  secondary  set  to  coincide  with  the  one 
perturbed  by  the  tides  on  the  primary  (so  r(t)  =  r*(t) ,  Msec  =  M*ec,  and  k  =  k*),  the 
above  expression  assumes  the  form  of 

5 

U(r)  &  -  \  GMsec  k2  cos  6.  (16b) 

4  r° 

Expressions  (14)  and  (16)  contain  the  longitudinal  lag 

e  =  ill  k*  At  =  2  (v*  -  8)  At  +  0(i2),  (17a) 


v*  and  6  being  the  true  anomaly  of  the  perturber  and  the  sidereal  angle  of  the  primary.  In  the 
special  case  (16b),  when  the  tide-generating  secondary  and  the  secondary  disturbed  by  the 
tide  on  the  primary  are  one  and  the  same  body,  the  asterisks  may  be  dropped: 


e  =  m  k  At  =  2(v  —  6)  At  +  0(i2). 


(17b) 


The  spin  rate  9  is  a  slow  variable,  in  that  it  may  be  assumed  constant  over  one  orbiting 
cycle.  The  true  anomaly  is  a  fast  variable.  For  a  nonvanishing  eccentricity,  v  too  is  a  fast 


5  From  r  =  a(  1  —  e2)/{\  +  e  cos  v)  and  dv/dM  =  (1  +  e  cos  v)2/(l  —  e 2)3/-  it  is  straightforward  that 


A r=r(t)  -  r(t  -  At) 


ae{  1  —  e2) 
(1+e  cos  v)2- 


sinv  Av+O  |e(Ar)2j 


ae  sin  v 
- -i  ■ n 

(1  -  e2)'/2 


A t+O  [e(n  AO2)  . 


Thus  our  replacement  of  r*(t  —  At)  with  r*(t)  =  r(t)  entails  a  relative  error  of  order  0(enAt).  In  expres¬ 
sions  (9),  (13),  (14),  and  (16),  the  absolute  error  will  be  of  the  same  order,  for  k  ^  X* .  However  for  X  =  k* 
the  absolute  error  will  become  0(enQ~^  At),  since  sin  €  is  of  the  same  order  as  the  inverse  quality  factor 
Q.  In  our  estimates  of  errors,  it  is  irrelevant  whether  we  define  Q  as  that  appropriate  to  the  principal  tidal 
mode  or  via  formula  (19)  below;  so  we  simply  use  the  generic  notation  Q.  As  we  explained  in  Efroimsky 
and  Williams  (2009),  after  averaging  over  one  revolution  of  a  nonresonant  secondary  about  the  primary,  the 
absolute  error  reduces  to  O  ^e^n^Q~ ^(Af)2^-  In  a  resonant  case,  though,  we  cannot  enjoy  this  reduction, 
because  in  this  case  the  averaging  procedure  looks  different  as  the  torque  changes  its  sign  over  a  period  of  the 
moon’s  revolution.  For  this  reason,  the  absolute  error  remains  0(enQ~l  At). 
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variable,  and  so  is  the  lag  6.  This  means  that  we  should  take  into  account  these  two  quantities’ 
variations  over  an  orbital  period. 

Had  we  not  changed  r3(f)  to  r3(f  —  At),  our  expression  (16b)  would  coincide,  up  to  an 
irrelevant  constant,  with  the  leading  term  of  the  appropriate  formula  from  MacDonald  (1964, 
Eq.  21).  To  appreciate  this  fact,  notice  that,  up  to  O  (z2),  the  absolute  value  of  the  longitudinal 
lag  e  is  the  double  of  the  geometrical  angle  A  =  \(6  —  u)  Af|  subtended  at  the  primary’s 
centre  between  the  directions  to  the  secondary  and  to  the  bulge,6  provided  At  is  postulated 
to  be  the  same  for  all  tidal  modes. 

The  analogy  between  the  MacDonald  theory  and  that  of  Darwin  and  Kaula  can  be  traced 
also  by  starting  from  the  series  (2).  Consider  the  case  of  the  tidally  perturbed  secondary 
coinciding  with  the  tide-raising  one,  so  X  =  X*,  and  all  the  orbital  variables  are  identical 
to  their  counterparts  with  an  asterisk.  It  will  then  be  easy  to  notice  that,  formally  (just  for¬ 
mally),  expression  (14)  mimics  the  principal  term  of  the  series  (2),  provided  in  this  term  the 
multiplier  G200  is  replaced  with  unity,  and  the  principal  phase  lag  e11QQ  =  2  (n  —  9)  Af 
is  replaced  with  the  longitudinal  lag  (17).  This  way,  within  the  MacDonald  formalism,  the 
longitudinal  lag  (17)  is  playing  the  role  of  an  instantaneous  phase  lag  associated  with  double 
the  instantaneous  synodic  frequency 

X  =  2\V  -  e\,  (18) 


which  is,  up  to  O  (i2),  the  double  of  the  angular  velocity  wherewith  the  point  located  under  the 
secondary  (with  the  same  latitude  and  longitude)  is  moving  over  the  surface  of  the  primary. 

To  extend  further  the  analogy  between  the  MacDonald  and  Darwin-Kaula  models,  one 
can  define  an  auxiliary  quantity 


"Q”  = 


(19) 


and  derive  from  (17)  and  (19)  that,  in  the  leading  order  of  e  and  z,  this  quantity  satisfies 


1 


X  A  t' 


(20a) 


A  popular  fallacy  would  then  be  to  interpret  (20a)  as  a  rheological  scaling  law  Q  ~  x“  with 
a  =  —  1 .  That  this  interpretation  is  generally  incorrect  follows  from  the  fact  that  the  quantity 
“<2”,  defined  through  (19),  is  not  obliged  to  coincide  with  the  quality  factor.7 


6  Be  mindful  that  the  double  of  the  geometrical  angle  is  not  equal  to  the  absolute  value  of  €9200  =  2(n  ~ 
0)  Af2200. 

7  Within  the  Darwin-Kaula  theory,  for  each  tidal  mode  o>impq,  we  introduce  the  phase  lag  as  eimpq  = 
^Impq  Atimpq .  Then  we  introduce  the  appropriate  quality  factor  Qimpq  =  Qiximpq )  =  2(KmP?l)  via 
the  expression  for  the  one-cycle  energy  loss: 


A Ecycle(Xtmpq )  — 


2  n  Epeak(xlinpq ) 

Q(X[mpq  ) 


Finally,  using  the  fact  that  Ximpq  =  \(,>tmpq  I  's  the  frequency  of  a  sinusoidal  load,  we  prove  that  the  afore 
introduced  and  Qimpq  are  interconnected  as  \/ Q  =  sin  e,  if  Epeak  denotes  the  maximal  energy,  or  in 

a  more  complex  way,  if  Epeak  stands  for  the  maximal  work  (Efroimsky  2012a, b). 

Within  the  MacDonald  method,  original  or  corrected  as  ( 15),  it  is  not  a  priori  clear  if  the  overall  peak  work 
(or  the  overall  peak  energy  stored)  and  the  overall  energy  loss  over  a  cycle  are  interconnected  via  the  auxiliary 
quantity  “Q”  in  exactly  the  same  manner  as  Epeak(Xj  )  and  A ECyCie(Ximpq)  are  interconnected  by 
QiXlmpq)  in  the  above  expression.  (Recall  that  the  total  cycle  of  the  tidal  load  is,  generally,  nonsinusoidal). 
Whenever  we  can  prove  that 
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To  sidestep  these  difficulties,  it  would  be  safer  to  write  the  constant-time-lag  rheological 
law,  for  small  lags,  simply  as 

|  6  |  =  x  At.  (20b) 

Historically,  MacDonald  (1964)  arrived  at  his  model  via  empirical  reasoning.  He  certainly 
realised  that  the  model  was  applicable  to  low  inclinations  only.  At  the  same  time,  this  author 
failed  to  notice  that  the  model  also  implied  the  frequency-independence  of  the  time  lag.  In 
fact,  this  frequency-independence,  (15),  is  necessary  to  derive  the  MacDonald  model  (16) 
from  the  generic  expression  (1)  for  the  tidal  amendment  to  the  primary’s  potential.  This  way, 
equality  (15),  is  a  priori  instilled  into  the  model.  In  other  words,  the  MacDonald  model 
of  tides  includes  in  itself  the  rheological  scaling  law  (20b)  with  a  constant  At.  Unaware 
of  this  circumstance,  MacDonald  (1964)  set  the  angular  lag  to  be  a  frequency-independent 
constant,  an  assertion  equivalent  to  the  time  lag  scaling  as  inverse  tidal  frequency.  However, 
as  we  just  saw  above,  a  consistent  derivation  of  the  MacDonald  tidal  model  requires  that  the 
time  lag  be  set  frequency-independent.8  Thus,  to  be  consistent,  the  MacDonald  method  must 
be  corrected  by  applying  (15)  or,  equivalently,  (20b).  In  greater  detail,  the  necessity  of  this 
amendment  is  considered  in  Efroimsky  and  Makarov  (2012). 

Even  after  the  model  is  combined  with  rheology  (15),  predictions  of  such  a  theory  are  of 
limited  use.  The  problem  is  that  this  rheology  is  radically  different  from  the  actual  behaviour 
of  solids.  As  a  result,  calculations  relying  on  (15)  render  implausibly  long  times  of  tidal  des¬ 
pinning — see,  for  example  the  discussion  and  references  in  Castillo-Rogez  et  al.  (201 1).  This 
tells  us  that  in  realistic  settings  the  MacDonald-style  approach  (14)  based  on  (15)  is  inferior, 
compared  to  the  Darwin  and  Kaula  method  (2),  which  may,  in  principle,  be  applied  to  any 
rheology.  Despite  this,  the  MacDonald  torque  remains  a  convenient  toy  model,  capable  of 
furnishing  results  which  are  qualitatively  acceptable  over  not  too  long  timescales  (Hut  1981; 
Dobrovolskis  2007). 


5  The  MacDonald  torque  and  the  ensuing  dynamical  model  by  Goldreich:  the  case 
of  an  oblate  body 

In  this  section,  we  shall  trace  how  the  MacDonald  theory  of  bodily  tides  yields  the  model 
of  near-resonance  spin  dynamics  by  Goldreich  (1966).  Then  we  shall  recall  an  oversight  in 
the  MacDonald  theory,  and  shall  demonstrate  that  correction  of  that  oversight  brings  a  minor 
alteration  into  Goldreich’s  model  of  spin  evolution. 


Footnote  7  continued 


p,( overall ) 

A  ^(overall )  peak 

AV  =  ^  ’ 

our  “<2”  can  be  spared  of  the  quotation  marks  and  can  be  called  instantaneous  quality  factor,  while  (20a) 
can  be  treated  as  a  reasonable  approximation  to  scaling  law  Q  ~  xa  with  a  =  —  1.  However,  this  should  be 
justified  in  each  particular  case,  not  taken  for  granted. 

8  In  application  to  a  non-resonant  setting,  a  constant- At  approach  was  taken  yet  by  Darwin  (1879).  Later, 
MacDonald  (1964)  abandoned  this  method  in  favour  of  a  constant-geometric  lag  calculation.  Soon  afterwards, 
though,  Singer  (1968)  advocated  for  reinstallment  of  the  constant- At  method.  Doing  so,  he  was  motivated  by 
an  apparent  paradox  in  MacDonald’s  treatment.  As  explained  by  Efroimsky  and  Williams  (2009),  the  paradox 
is  nonexistent.  Nonetheless,  the  work  by  Singer  (1968)  was  fruitful  at  the  time,  as  it  renewed  the  interest  in  the 
constant-At  treatment.  The  full  might  of  the  model  was  revealed  by  Mignard  (1979,  1980,  1981)  who  used  it 
to  develop  closed  expressions  for  the  tidal  force  and  torque. 
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The  goal  of  this  section  is  limited  to  consideration  of  the  tidal  torque  solely.  So  we  assume 
that  the  body  is  oblate,  and  the  triaxiality-caused  torque  does  not  show  up.  Appropriate 
for  spinning  gaseous  objects,  this  treatment  indicates  that  the  tidal  torque  stays  finite  for 
synchronous  rotation  and  vanishes  at  an  angular  velocity  slightly  faster  than  synchronous. 

In  Sect.  6  below  we  shall  address  the  more  general  setting  appropriate  to  telluric  bodies, 
with  triaxiality  included.  While  some  authors  (e.g.,  Heller  et  al.  201 1)  ignore  the  triaxiality- 
caused  torque  in  their  treatment  of  solid  planets,  it  turns  out  that  inclusion  of  the  permanent- 
figure  torque  renders  important  physical  consequences  and  changes  the  picture  completely, 
as  will  be  seen  in  Sect.  6. 

5.1  The  MacDonald  torque 

We  shall  restrict  ourselves  to  the  case  of  the  tidally  perturbed  secondary  coinciding  with  the 
tide-raising  one,  so  Msec  =  M*ec,  and  all  the  orbital  variables  are  identical  to  their  counterparts 
with  an  asterisk. 

Differentiating  (14)  with  respect  to  X,  and  then  setting  X  =  X*,  we  obtain  the  following 
expression  for  the  polar  component  of  the  torque,  for  low  i: 

5 

Tz  =  l  GM2ec  k2  — - * sin  (2(v  -  6)  Af)  +  Od2/Q) 

1  r  (t)  r*  (t-  At) 

3  R 

=  -  GM;ec  k2  —  sin  (2(v  -  0)  At)  +  0(i2/Q )  +  0(enAt/Q ),  (21) 

2  r 

where  the  error  0(enAt/Q )  emerges  when  we  identify  the  lagging  distance  r*{t  —  At) 
with  r*(t)  =  r(t),  as  explained  in  footnote  5  in  the  preceding  section.  Referring  the  Reader 
to  Efroimsky  and  Williams  (2009)  for  this  and  other  technicalities,  we  would  mention  that 
(21)  is  equivalent  to  the  Darwin  torque  only  under  the  condition  that  the  rheological  model 
(15-20a)  is  accepted.  For  potentials,  employment  of  model  (15-20a)  enables  one  to  wrap  up 
the  infinite  series  (2)  into  the  elegant  finite  form  (9).  For  torques  (truncated  to  l  =  2  only), 
similar  wrapping  of  the  appropriate  series  is  available  within  the  said  model. 

In  the  preceding  section,  we  explained  the  geometric  meaning  of  the  longitudinal  lag  (17): 
its  absolute  value  is  the  double  of  the  geometric  angle  separating  the  directions  to  the  bulge 
and  the  secondary  as  seen  from  the  primary’s  centre. 

If  we  define  a  quantity  “Q”  via  (20a),  the  MacDonald  torque  will  look: 

5 

Tz=3-  GM2c  k2— - * -  sin 6  +  0(i2 1 Q) 

z  r  (t)  r  *  (t  —  At) 

5 

=  |  GM2c  k2  A.  _i_  Sgn(u  -6)+  0(i2/ Q)  +  0(enAt/Q).  (22) 

For  a  nonzero  eccentricity,  the  quantity  “Q”  should  not  be  interpreted  as  an  instantaneous 
quality  factor,  because  it  is  not  guaranteed  to  interconnect  the  peak  work  or  peak  energy  and 
the  one-cycle  energy  loss  in  a  manner  appropriate  to  a  quality  factor — see  footnote  7  in  the 
previous  section.  Therefore  a  more  reasonable  and  practical  way  of  writing  the  MacDonald 
torque  (21)  would  be  through  using  (20b)  or  (17): 
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3  R~ 

Tz  =  -  GMjec  k2—  At  2  (v  -  9)  +  0(i2/Q )  +  0{enAt/Q).  (23) 

-  r 

As  we  saw  above,  the  MacDonald  model  is  self-consistent  only  for  At  being  a  frequency- 
independent  constant.  However  the  factor  v  —  9  showing  up  in  (23)  varies  over  a  cycle,  for 
which  reason  the  torque  needs  averaging.  This  averaging  is  carried  out  in  Appendix  B.  In 
the  vicinity  of  the  1:1  resonance  (for  9  close  to  n),  the  sign  of  v  —  9  changes  twice  over  a 
cycle,  which  makes  the  averaging  procedure  nontrivial.  This  situation  is  to  be  addressed  in 
Sects.  5.2  and  5.3. 

5.2  Goldreich  (1966):  treatment  based  on  the  MacDonald  torque 


In  this  subsection,  we  shall  briefly  recall  a  presently  conventional  method  pioneered  almost 
half  a  century  ago  by  Goldreich  (1966).  Based  on  the  MacDonald  tide  theory,  this  method 
has  inherited  both  its  simplicity  and  its  flaws. 

The  1:1  resonance  takes  place  when  the  spin  rate  9  of  the  primary  (the  satellite)  is  equal 
to  the  mean  motion  n  wherewith  the  secondary  body  (the  planet)  is  apparently  orbiting  the 
primary.  The  formula  (22)  for  the  MacDonald  torque  contains  not  the  difference  9  —  n  but 
the  difference  9  —  v,  for  which  reason  the  expression  under  the  integral  may  twice  change 
its  sign  in  the  course  of  one  revolution. 

With  aid  of  the  formula 


v  =  M  +  2e  sin.M  +  —  e2  sin2AI  +  0(e 3)  (24) 

and  under  the  assumption  that  e  «  n,  we  shape  the  difference  of  our  concern  into  the  form 
of 


9  —  v  =  (9  —  n)  —  2ne  ^cosAf  +  cos2.M^  +  0(e3) 
=  fj  —  2  n  e  ^cosAJ  +  ^ecos2A1^  +  0(e3) 


or,  equivalently, 

9  -  v 


2/i  e 


cos  A!  +  -  e  cos2A1  —  — - —  +  0(e 2) 
4  2  n  e 


(25) 


(26) 


Here 


rj  =  9  —  M  —  co  —  Si  (27) 

is  a  slowly  changing  quantity,  whose  time-derivative  r)  =  9  —  n  becomes  nil  when  the 
system  goes  through  the  resonance.9  To  impart  the  words  “slowly  changing”  with  a  definite 
meaning,  we  assert  that  i)/n  is  of  order  e 2 — a  claim  to  be  justified  a  posteriori. 

Expression  (25)  changes  its  sign  at  the  points  where 

cosTVI  =  —  -  e  cos  2 M  +  — - —  +  0(e2).  (28) 

4  2  n  e 


9  Goldreich  (1966)  defined  this  quantity  simply  as  =  9  —  A4,  because  he  reckoned  8  from  a  fixed  perihe¬ 
lion  direction.  We  however  reference  our  6  from  a  direction  fixed  in  space.  (It  is  the  same  direction  wherefrom 
the  node  is  reckoned).  This  convention  originates  from  our  definition  of  9  as  the  sidereal  angle — it  is  in  this 
capacity  that  9  was  introduced  back  in  Eqs.  (3-4).  As  we  are  not  considering  the  nodal  or  apsidal  precession, 
our  subsequent  formulae  containing  9  will  be  equivalent  to  those  ensuing  from  Goldreich's  definition  of  9. 
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As  all  the  terms  on  the  right-hand  side  of  (28)  are  of  order  e  at  most,  so  must  be  the  term  on 
the  left-hand  side.  Hence  condition  (28)  is  obeyed  in  the  two  points  whose  mean  anomaly 
(and  therefore  also  true  anomaly)  is  close  to  ±  jr/2: 

V  =  ±  (I  “  S)’ 

S  being  of  order  e.  From  (29)  and  (24)  we  obtain: 

sin  8  =  cos  v  =  cos(Af  +  2  e  sin  M  +  0(e2)) 

=  cosA4  —  2e  sin2 *  TV)  +  0(e2) 
which,  in  combination  with  (28),  entails: 

sin  S  = - e  cos  2 M  +  — —  —  2  e  sin2  M  +  O (e2) 

4  2  n  e 

=  ~  e  (5-  -  \  sin2Af)  +  0(e2).  (31) 

2  ne  \  4  2  ) 


(29) 


(30) 


Insertion  of  (29)  into  (24)  also  yields  ±  sinAI  =  cos  5+  0(e).  Recalling  that  S  is  of  order 
e,  we  obtain:  sin2  M  =  1  —  sin2  8  +  0(e)  =  1  +  0(e).  This  enables  us  to  rewrite  (31) 
as 

-  \  e  +  0(e2).  (32) 

2  ne  4 

Before  finding  the  rate  f]  =  9  —  n  at  which  the  resonance  is  traversed,  let  us  enquire  if 
perhaps  it  could  be  simply  put  nil,  the  satellite  being  permanently  kept  in  the  resonance.  The 
answer  is  negative,  because  for  a  vanishing  6  —  n  the  difference  9  —  v  emerging  in  (25) 
becomes  a  varying  quantity  of  an  alternating  sign,  and  so  becomes  the  torque.  On  general 
grounds,  one  should  not  expect  that  the  average  torque  becomes  nil  for  a  vanishing  ?),  though 
it  may  vanish  for  some  finite  value  of  i). 

To  find  this  value  of  r),  many  authors  (Goldreich  1966,  Eq.  15;  Kaula  1968,  Eq.  4.5.29; 
Murray  and  Dermott  1999,  Eq.  5.11)  simply  integrated  the  MacDonald  torque  (22),  assuming 
the  quality  factor  Q  constant,  and  thus  keeping  it  outside  the  integral: 


(Tz) 


( Goldreich ) 


3  GM~k2R 


4  tt  a2  Q  (l  —  e2)1,/2  •'  r 


2tt  4 

r  r 

~r  sgn(0  ~  v)  dv 


3  GM  k2  R 


4  7r  a-  Q  (i  _  e2) 


1/2 


it/2—8  4 

R 


R 


— r  dv  —  /  — r  dv 


3  GMsec  k2  R5 
2  jr  a6  Q 


jt/2—8  ji 


tt/2—8 


ti/2-8 


(1  +  4e  cos  v)  dv  +  0(e~) 


3  GM„r  ki  R 5  , 

=  - “6C  2 - (4e  cos  5  -  5)  +  0(e~). 

n  a°  Q 

In  a  trapping  situation,  the  average  torque  vanishes.  So  (34)  entails: 

S  =  4e  cos  S  =  4  e  +  0(e2). 


(33) 


(34) 


(35) 
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By  combining  the  latter  with  (32),  the  afore  quoted  authors  arrived  at 

19  2 

n stall  =  ~ne  ,  (36) 

an  expression  repeated  later  in  papers  and  textbooks  (e.g.,  Eq.  4.5.37  in  Kaula  1968,  or 
Eq.  5.14  in  Murray  and  Dermott  1999).  This  result  however  needs  to  be  corrected,  because 
the  MacDonald  approach  is  incompatible  with  the  frequency-independence  of  Q  assumed  in 
(34). 

Be  mindful  that  (36),  as  well  as  its  corrected  version  (44)  to  be  derived  below,  indicate 
that  r]  is  of  order  e2.  This  justifies  our  assertion  made  in  the  paragraph  after  formula  (27). 

5.3  The  corrected  MacDonald  model 


To  impart  the  MacDonald  treatment  consistently,  one  has  to  calculate  the  averaged  torque, 
with  the  frequency-dependence  of  Q  taken  into  consideration.  As  we  saw  in  Sect.  4,  the 
MacDonald  approach  fixes  this  dependence  in  a  manner  that  can,  with  some  reservations,  be 
approximated  with 


or,  in  more  general  notations, 


1 

X  A  t 


(37) 


Q  =  £a  x“,  with  a  =  -  1, 


(38) 


where  the  double  instantaneous  synodic  frequency  is  given  by  (18). 

The  form  (38)  of  the  rheological  law  is  more  convenient,  as  it  leaves  us  an  opportunity 
to  consider  values  of  a  different  from  —  1 .  For  any  value  of  a,  the  constant  £  is  an  integral 
rheological  parameter,  which  has  the  dimension  of  time,  and  whose  physical  meaning  is 
discussed  in  Efroimsky  and  Lainey  (2007).  It  can  be  demonstrated  that  in  the  special  case 
of  a  =  —  1  the  parameter  £  coincides  with  the  time  lag  At.  For  actual  terrestrial  bodies,  a 
is  different  from  —1,  and  the  integral  rheological  parameter  £  is  related  to  the  time  lag  in  a 
more  complicated  manner  (Ibid.). 

As  demonstrated  in  Appendix  B,  insertion  of  (38)  and  (18)  into  (22),  for  a  =  —  1 ,  or  equiv¬ 
alently,  direct  employment  of  (23),  entails  the  following  expression  for  the  orbit-averaged 
torque  acting  on  a  secondary: 


(Tz)  =  -  Z  [  6  A(e)  -  n  N(e)  ]  +  0(i2/Q)  +  0(Q~3)  +  0(enAt/Q),  (39) 


where 


A(e)  —  I  1  +  3  e  +  —  £ 


)  (l-eT"'2 


15  ,  105  ,  , 

=  1  +  —  e2  +  —  e4  +  0(e6)  (40) 


and 


15  .  45 


-vw  =  (1  +  r  +  i'  +  tA 

+  0(e6), 

while  the  factor  Z  is  given  by10 


’)  0-eT 


,  27  2  573  4 

=  1  +  —  ^  +  —  e 


(41) 


10  Recall  that,  for  a  =  —  1,  the  rheological  parameter  E  is  simply  the  time  lag  At. 
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z  =  3  GM]ec  k2  £  =  3 n2M*ec  k2  At  R '5_  =  3 n  M*c  k2  R ,5_  n_ 

R  a, ^  ( Mprim  H"  Msec)  Q  ( Mpr{m  +  Msec)  Cl ^  X 

M prim  and  Msec  being  the  masses  of  the  primary  and  the  secondary. 

Be  mindful  that  the  right-hand  side  of  (42)  contains  a  multiplier  y  =  - — - which 

2\6  -v  | 

is  missing  in  the  despinning  formula  employed  by  Correia  and  Laskar  (2004, 2009).  This  hap¬ 
pened  because  in  Ibid,  the  quality  factor  was  introduced  as  1  /(n  At)andnotas  l/(x  Af) — see 
the  line  after  formula  (9)  in  Correia  and  Laskar  (2009).  In  reality,  the  quality  factor  Q  must, 
of  course,  be  a  function  of  the  forcing  frequency  /  (which  happens  to  coincide  with  the  mean 
motion  n  in  the  3:2  and  1:2  resonances  but  differs  from  n  outside  these). 

The  quality  factor  Q  being  (within  this  model)  inversely  proportional  to  x ,  the  presence 
of  the  y  factor  in  (42)  makes  the  overall  factor  Z  a  frequency-independent  constant. 

A 

Rewriting  (39)  as 


(Tz) 


=  -Z 


H 

f/Q 

["( 


=  -z 

+  Q{i2/Q)  +  0(enAt/Q) 


15 


1  +  /f-  e~  +  —  e4  |  - 


105 


)  -  6~2  (■ + w  el)\ 


+  0(e°) 


+0(i1/Q)+0(enAt/Q), 


(43) 


we  see  that  it  vanishes  when  the  rate  of  change  of  r]  =  8  —  M  —  u>  —  £2  accepts  the  value 


fj  t  ,,  =  6  n  e1  +  -  n  e4  +  0{eb)  +  0{i~/Q)  +  0(enAt/Q). 


(44) 


For  r)  larger  or  smaller  than  f]stau ,  the  average  torque  (43)  is  nonzero  and  impels  i)  to  evolve 
towards  the  stall  value  (44). 

On  the  right-hand  side  of  (44),  the  leading-order  term  contains  a  numerical  factor  of  6,  as 
different  from  the  factor  19/2  showing  up  in  (36).  The  necessity  to  change  19/2  to  6  in  the  e 2 
term  was  pointed  out  by  Rodriguez  et  al.  (2008,  Eq.  2.4),  who  had  arrived  at  this  conclusion 
through  some  different  considerations  [which,  too,  were  based  on  the  frequency-dependence 
(38)].  The  same  result  can  be  obtained  within  the  Darwin-Kaula  approach,  provided  the 
scaling  law  (38)  is  employed  (Efroimsky  2012a,b).n 

Although  correction  of  the  oversight  in  the  MacDonald  torque  renders  a  number  different 
from  the  one  furnished  by  Goldreich’s  development  (6  instead  of  19/2),  qualitatively  the 
principal  conclusion  by  Goldreich  (1966)  remains  unchanged:  when  an  oblate  body’s  spin 
is  evolving  toward  the  resonance,  vanishing  of  the  average  tidal  torque  entails  spin  slightly 
faster  than  resonant,  a  so-called  pseudosynchronous  rotation.  It  however  should  be  strongly 
emphasised  that  the  possibility  of  pseudosynchronous  rotation  hinges  upon  the  dissipation 
model  employed.  Makarov  and  Efroimsky  (2012)  have  found  that  a  more  realistic  tidal  dis¬ 
sipation  model  than  the  corrected  MacDonald  torque  makes  pseudosynchronous  rotation 
impossible. 


11  Implicitly,  this  result  is  present  also  in  Correia  et  al.  (2011,  eqn  20),  Laskar  and  Correia  (2004),  eqn  9,  and 
in  Hut  (1981).  The  earliest  implicit  occurrence  of  this  result  was  in  Goldreich  and  Peale  (1966,  eqn  24). 
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6  Evolution  of  rotation  near  the  1:1  resonance:  the  case  of  a  triaxial  body 

In  distinction  from  a  gaseous  or  liquid  body,  a  solid  body  would  be  expected  to  have  a  perma¬ 
nent  figure  in  addition  to  the  tidal  distortion  discussed  so  far.  Goldreich  (1966)  demonstrated 
that  this  permanent  figure  plays  an  important  role  in  determining  if  a  primary,  whose  rotation 
is  being  slowed  down  or  sped  up  by  the  tidal  torque  caused  by  a  secondary,  can  be  captured 
into  the  synchronous  rotation  state.  This  section  follows  Goldreich’s  derivation,  but  substi¬ 
tutes  the  MacDonald  tidal  torque  with  its  corrected  version,  and  also  uses  a  more  general 
mass  expression. 

6.1  Rotating  primary  subject  to  a  triaxiality-caused  torque:  the  equation  of  motion 
and  the  first  integral 

Consider  a  triaxial  primary  body,  which  has  its  principal  moments  of  inertia  ordered  as 
A  <  B  <  C ,  and  which  is  rotating  about  the  maximal-inertia  axis  associated  with  moment 
C.  About  the  primary,  a  secondary  body  describes  a  near-equatorial  orbit  (so  its  inclination  on 
the  primary’s  equator,  i,  may  be  neglected).  The  secondary  exerts  on  the  primary  two  torques. 
One  being  tidal,  the  other  is  triaxiality-caused,  i.e.,  generated  by  the  existence  of  the  perma¬ 
nent  figure  of  the  primary.  Its  component  acting  on  the  primary  about  the  maximal-inertia 
axis  is 

3  GMsec 

Ttria*  =  ~  (B  ~  A)  sin  2X,  (45) 

Z  I 

the  longitude  X  being  furnished  by  formula  (10).  In  that  formula,  the  sidereal  angle  9  is 
reckoned  from  a  reference  direction  in  space  to  the  principal  axis  associated  with  moment  A. 
The  acceleration  of  the  sidereal  angle  then  obeys 

..  3  GMsec 

C  9  -  -  (B  -  A)  - ^  sin  2X  =  0.  (46) 

2  ri 

In  neglect  of  the  nodal  and  apsidal  precession,  as  well  as  of  Ad0,  definition  (27)  yields: 12 

■h  =  6  -  n.  (47) 

Ignoring  changes  in  the  mean  motion,13  we  write: 

9  =  ij,  (48) 

so  the  first  term  in  (46)  becomes  simply  Cij. 


12  The  caveat  about  three  neglected  items  implies  that  our  n  is  the  osculating  mean  motion  n(t)  =  y/  p/a (t)3, 
and  that  we  extend  this  definition  to  perturbed  settings.  The  so-defined  mean  anomaly  evolves  in  time  as 
Ai  =  _Mo(f)  +  f,  n(t)dt,  whence  A4  =  Aig  +  n(t). 

The  said  caveat  becomes  redundant  when  n  is  defined  as  the  apparent  mean  motion,  i.e.,  either  as  the  mean- 
anomaly  rate  djvl/dt  or  as  the  mean-longitude  rate  dL/dt  =  dO./dt  +  dw/dt  +  dAd/dt  (Williams  ct  al. 
2001).  While  the  first-order  perturbations  of  a(t)  and  of  the  osculating  mean  motion  •/ p./ait )-  include  no 
secular  terms,  such  terms  are  often  contained  in  the  epoch  terms  tl ,  do,  and  -Mo-  This  produces  the  differ¬ 
ence  between  the  apparent  mean  motion  defined  as  dL/dt  (or  as  dAd/dt)  and  the  osculating  mean  motion 

V/r/a(f)3. 

Our  neglect  of  evolution  of  the  mean  motion  is  acceptable,  because  orbital  acceleration  h  is  normally  much 
smaller  than  the  spin  acceleration  6.  Indeed,  for  torques  arising  from  the  primary,  h/6  is  of  the  same  order  as 
the  ratio  of  C  to  the  orbital  moment  of  inertia. 
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To  process  the  second  term  in  (46),  we  would  compare  (10)  with  (27): 
a  =  —  (9  4~  T"  co  +  v  O  (i2')  =  (  —  0  +  4-  co  4~  Al)  4~  (v  —  Af)  4~  O  (/“) 

=  -  r>  +  (v  -  M)  +  0(i 2)  (49) 

In  neglect  of  the  inclination,  the  following  approximation  is  acceptable  in  the  vicinity  of  the 
1:1  resonance:14 

(r-3sin2A)  =  —  G2oo(e)  a~3  sir\2ri,  (50) 

where  (•  •  • )  signifies  orbital  averaging,  while  the  eccentricity  function  can  be  approximated 
with 

G2 oo(e)  =  1  -  -  e2  +  0(e4).  (51) 

This  way,  omitting  0(i2)  in  (49)  and  substituting  sin  2k  with  its  average  in  (46),  we  transform 
the  second  term  in  (46)  to: 

3  GMsec 

+  (B  -  A)  - ^  G2 oo(e)  sm2ij.  (52) 

2 

While  Goldreich  (1966)  assumed  that  the  mass  of  the  secondary  is  much  larger  than  the  mass 
of  the  primary,  we  do  not  impose  this  restriction.  Combining  Kepler’s  third  law,  G{Msec  4- 
Mprim)/a3  =  n2,  with  formulae  (46),  (48),  and  (52),  we  finally  arrive  at 

Cij  +  l  (B  -  A)  Msec  - n2  G200  (e)  sin  217  =  0,  (53) 

—  Msec  4“  Mpr  im 

an  equation  describing  the  evolution  of  the  rotation  angle  rj.  As  pointed  out  by  Goldreich 
(1966),  this  equation  is  equivalent  to  the  one  describing  a  simple  pendulum.  Indeed,  in  terms 
of  B  =  2 n,  Eq.  (53)  becomes  B  +  y2  sin  B  =  0,  with  a  constant  positive  y2  and 
with  /3  playing  the  role  of  the  pendulum  angle.15 

It  follows  directly  from  (53)  that  the  spin  acceleration  ij  vanishes  if  r]  assumes  the  values 
of  0  or  7r.  As  can  be  seen  from  the  pendulum  analogy,  the  initial  conditions  (77 ,  f])  t=t  = 
(0  ,  0),  as  well  as  the  conditions  (??,)))  t=t  =  {n  ,  0),  correspond  to  the  situation  where 
the  pendulum  comes  to  a  stall  in  the  lower  point  (so  ji  =  0  when  fi  =  0  or  2tt).  Under  such 
initial  conditions,  rj  stays  0  or  n  all  the  time,  which  implies  a  uniform  synchronous  rotation 
of  the  secondary  about  the  primary.16 

Multiplication  of  Eq.  (53)  by  r],  with  subsequent  integration  over  time  t,  gives  the  first 
integral  of  motion, 

\  C  >72  ~-AB  -  A)  Msec  - n2  G2 00(e)  cos  2ri  =  E,  (54) 

2  4  MSec  4“  Mprim 

whose  value  depends  on  the  initial  conditions. 

14  Pioneered  by  Goldreich  and  Peale  (1966),  the  approximation  is  explained  in  more  detail  by  Murray  and 
Dermott  (1999)  whose  treatment  omits  terms  of  the  order  of  e ^  and  higher  (see  formulae  5.59-5.60  in  Ibid.). 
To  justify  the  omission,  recall  that  before  averaging  the  expansion  of  r-3  sin  2  X  includes  a  series  of  terms 
with  different  arguments  of  the  form  sin(2?7  +  qA4),  where  q  is  an  integer.  Since  we  are  interested  in  dynamics 
in  the  vicinity  of  the  1:1  resonance,  where  q  is  a  slow  variable,  then  only  the  sin(2^)  term  remains  after  the 
orbital  averaging.  The  averaged-out  terms  are  of  the  order  of  e  and  higher  powers  including  the  e 3  terms. 

15  In  Sect.  6.5  below,  we  write  down  the  expression  for  the  libration  frequency  Xnb—max  and  also  explain 
the  reason  why  we  equip  it  with  such  a  subscript — see  formulae  (68)  and  (69). 

16  According  to  (53),  the  spin  acceleration  ij  vanishes  also  for  q  =  ±n/2,  which  is  the  upper  point  of  the 
pendulum. 
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The  period  of  the  variable  ij  is  given  by  a  quadrupled  integral  over  a  quarter-cycle  of  rj: 


max 


r)= 0 


For  circulation,  use  rjmax  =  n 1 2.  In  the  case  of  libration  about  ij  =  0,  the  value  of  rjmax  is 
less  than  n / 2  and  can  be  expressed  via  E  by  setting  rj  =  0  in  Eq.  (54).  To  get  an  explicit 
expression  of  P,  one  should  first  express  rj  via  E  using  (54)  and  taking  the  positive  root,  and 
then  should  plug  the  so-obtained  expression  for  i)  into  (55).  Also  be  mindful  that  libration 
about  7T  can  be  converted  to  the  same  integral  by  adding  or  subtracting  j r  from  the  variable 
of  integration. 

In  what  follows,  averaging  over  the  period  P  will  be  denoted  by  (•  •  • ) p,  the  subscript 
serving  to  distinguish  the  operation  from  averaging  over  the  orbit  employed  in  Sect.  5  and  in 
formula  (50). 


6.2  Parallels  with  pendulum:  the  auxiliary  quantity  W 


Goldreich  (1966)  noted  the  similarity  of  the  differential  equations  (53-54)  to  the  classical 
pendulum  problem.  If  the  initial  conditions  place  the  body  outside  of  the  1 : 1  spin-orbit  reso¬ 
nance,  then  i]  circulates  with  a  forced  oscillation  in  rotation  that  depends  on  the  mean  value 
{rj)  p  of  the  fi  frequency.  If  however  the  initial  conditions  place  the  body’s  spin  within  a 
sufficiently  close  proximity  of  the  1 : 1  resonance  (the  “resonance  region”),  then  rj  will  librate 
about  0  or  jt  .  This  will  be  a  free  physical  libration  with  an  amplitude  smaller  than  7r/2.  Inside 
the  resonance  region,  the  initial  conditions  establish  the  free-libration  amplitude  and  phase. 
However  the  restricted  nature  of  the  motion  will  keep  the  mean  value  of  tj  constant:  it  will  be 
either  0  or  jt  .  For  the  same  reason,  (rj) p  will  remain  zero  in  the  libration  regime.  In  contrast 
to  this,  outside  of  the  resonance  region  the  mean  rate  of  circulation  (rj)  p  will  possess  a  value 
determined  by  the  initial  conditions  on  ij  and  rj. 

Employing  (54)  to  express  rj  via  E,  and  plugging  this  expression  into  (55),  Goldreich 
(1966)  demonstrated  that  the  period  P  can  be  expressed  via  E,  inside  or  outside  the  reso¬ 
nance  region,  by  a  complete  elliptic  integral  of  the  first  kind.  This  is  natural,  as  the  expression 
of  ij  through  E  involves  a  square  root. 

Just  as  in  the  pendulum  case,  there  exists  a  critical  E  for  which  the  period  diverges.  This 
is  the  boundary  E  =  Ef,  separating  circulation  from  libration.  To  locate  the  boundary,  one 
should  set  simultaneously  rj  =  jt/2  and  rj  =  0  in  formula  (54): 


A) 


Msec 

Msec  T"  Mpyim 


G2oo(e)  «2 • 


(56) 


Mathematically,  the  emergence  of  a  logarithmic  singularity  in  P  at  E  =  E j,  can  be  observed 
from  formulae  (9—10)  in  Ibid.  Physically,  this  situation  resembles  the  slowing-down  of  a  pen¬ 
dulum  at  the  circulation/libration  border.  In  our  problem,  though,  this  division  corresponds 
to  rj  =  ±7r/2. 

By  setting  simultaneously  rj  =  0  and  rj  =  0  in  (54),  we  obtain  the  minimal  value  that  the 
integral  E  can  assume: 


min  E  =  —  E/,. 


(57) 
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The  values  E  >  E\,  correspond  to  circulation  in  either  direction,  with  oscillations  of  the 
rotation  rate.  The  values  falling  within  the  interval  Eh  >  E  >  —  Eb  give  libration.  The 
minimum  value  E  =  —  E\,  corresponds  to  synchronous  rotation  without  libration. 

The  quantity  r\  being  slowly  varying  (compared  to  the  mean  motion  /?),  the  period  P, 
with  which  ;;  is  changing,  naturally  turns  out  to  be  much  longer  than  the  orbital  period,  both 
outside  and  inside  the  1:1  resonance: 

2n 

P  »  — ■  (58) 

n 

This  can  be  appreciated  from  the  evident  equalities 

outside  of  the  1 : 1  resonance,  (59a) 

inside  the  1:1  resonance.  (59b) 


P  = 


P  = 


2  7T 

in )P  I 

2  7T 

Xlib 


In  (59b),  Xiib  denotes  the  libration  frequency  which  too  is  much  smaller  than  the  mean 
motion  n,  as  we  shall  see  shortly. 

A  useful  quantity  defined  by  Goldreich  (1966)  was 


W  =  P(r )2)P 


V—^max 


rj=0 


(60) 


Similarly  to  (55),  r)max  =  tt /2  for  circulation,  and  i]max  <  jt/2  for  libration.  Just  as  P,  so  W 
can  be  expressed  through  E  by  complete  elliptic  integrals.  For  circulating  (f/) p,  Goldreich 
described  P  and  W  as  applying  to  one  oscillation,  but  they  describe  one  circulation  of  t]  with 
two  oscillations.  For  librations,  P  and  W  describe  one  libration  cycle. 

One  more  quantity  of  use  in  this  problem  will  be  the  mean  square  variation  of  rj  about 
{fi) P,  given  by  (rj2) p  —  ( ?) )2 .  Be  mindful  that  the  notation  {•••)/>  is  employed  to  indicate 
averaging  over  a  circulation  cycle  as  well  as  that  over  a  libration  cycle. 

When  the  values  of  the  mean  motion  n,  the  mass  factor  Msec/(Msec  +  Mprim),  and  the 
ratio  (B  —  A)/ C  are  given,  and  the  initial  conditions  on  7/  and  r)  are  set,  Eq.  (54)  furnishes 
the  value  of  E/C.  The  comparison  of  this  value  with  Eb/C  distinguishes  circulation  from 
libration.  With  the  values  of  (B  —  A)/C  and  E / C  known,  the  complete  elliptic  integrals  can 
be  evaluated  and  the  values  of  the  quantities  P  ,  W  ,  (?))  p ,  and  ( f} 2)  p  can  be  found.  It  should 
also  be  possible  to  reverse  this  procedure.  For  example,  the  knowledge  of  the  circulating 
value  of  (r/) p,  along  with  the  knowledge  of  n,  the  mass  factor  Msec/{Msec  +  Mpr!m),  and 
the  ratio  (B  —  A)/ C,  should  allow  E/C  ,  W,  and  (j)2) p  to  be  computed. 


6.3  The  tidal  torque  and  the  libration  bias 


To  account  for  tidal  dissipation,  Goldreich  (1966)  added  the  averaged  (over  an  orbital  period) 
tidal  torque  (34)  to  the  right-hand  side  of  Eq.  (53).  We  however  shall  employ  the  corrected 
average  torque  (43)  instead  of  (34).  This  will  result  in 


Crj  +  \{B  -  A)  — — n2  G2 oo(e)  sin2?j 

2  MSec  Mprim 

=  -,[,(1+^2)-6we2(i+f  e2) 

+  0(enAt /Q), 


+  0{e 6)  +  0(i2/g) 


(61) 
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where  the  coefficient  Z  depends  via  (42)  on  the  orbital  variables,  the  tidal  parameters,  and 
the  masses. 

As  we  already  mentioned  at  the  beginning  of  Sect.  5,  ignoring  the  permanent-figure  term 
would  lead  one  to  the  conclusions  that  the  tidal  torque  is  finite  for  synchronous  rotation  and 
that  it  vanishes  for  a  spin  rate  slightly  higher  than  synchronous.  However  inclusion  of  the 
permanent-figure  term  into  the  picture  alters  the  results  radically.  The  constant  term  on  the 
right-hand  side  of  Eq.  (61)  causes  the  mean  value  of  r)  to  be  slightly  larger  than  0  or  it .  For 
small  librations,  this  bias  is 


+  0(i2/Q )  +  0(enAt/Q) 


(62a) 


+  0(i2/Q)  +  0{enAt/Q ) 


(62b) 


with  Z  calculated  through  (42).  The  value  of  the  time  lag  At  entering  the  expression  (42) 
for  Z  should  be  set  appropriate  to  the  orbital  frequency  (mean  motion). 

As  evident  from  (62),  the  bias  comes  into  being  due  to  the  eccentric  shape  of  the  orbit. 
Bodies  with  permanent  figures  can  achieve  synchronous  rotation  because  the  bias  in  i]  gives 
birth  to  a  permanent-figure  torque  that  balances  the  constant  dissipative  torque. 

At  this  point,  an  important  caveat  will  be  in  order.  As  we  explained  in  Sect.  4,  the  MacDon¬ 
ald  model  of  tides  becomes  self-consistent  only  when  the  time  delay  At  emerging  in  (42)  is 
set  frequency-independent.  This  circumstance  limits  the  precision  of  the  MacDonald  torque 
and  of  the  Goldreich  dynamical  model  based  thereon,  whenever  the  model  gets  employed  to 
determine  the  timescales  of  spin  evolution  (Efroimsky  and  Lainey  2007).  Nevertheless,  for 
very  slow  evolution  the  model  may  be  employed  for  obtaining  qualitative  estimates.  In  this 
case,  At  should  be  treated  as  a  parameter  that  itself  depends  upon  the  forcing  frequency  in  the 
material.  For  rotation  outside  the  1 : 1  spin-orbit  lock,  it  would  be  a  tolerable  approximation 
to  use  At  appropriate  to  the  principal  tidal  frequency  X2200  or  to  double  instantaneous 
synodic  frequency  (18).  However  inside  the  1:1  resonance.  At  would  correspond  to  the 
libration  frequency  Xub  which  may  be  very  different  from  the  usual  tidal  frequencies  for 
nonsynchronous  rotation.  This  circumstance  may  change  the  value  of  At  and  therefore  of  Z 
noticeably. 

6.4  Tidal  dissipation  and  the  point  of  stall 

Be  mindful  that  on  the  right-hand  side  of  Eq.  (61)  we  have  the  tidal  torque  averaged  over 
the  orbit  (see  Appendix  B  for  details).  Similarly,  the  second  term  on  the  left-hand  side  of 
(61)  is  the  permanent-figure  torque  averaged  over  the  orbit — see  expression  (52). 17  There¬ 
fore  Eq.  (61)  renders  us  the  behaviour  of  r]  over  times  longer  than  the  orbital  period.  This 
is  acceptable,  because  the  orbital  period  is  much  shorter  than  the  timescales  of  our  interest. 
Specifically,  it  is  much  shorter  than  P ,  see  Eq.  (58). 

Without  dissipation,  E  is  conserved  during  oscillations  of  77,  while  in  the  presence  of  weak 
dissipation  E  changes  slowly.  Indeed,  multiplying  both  sides  of  (6 1 )  by  77  and  making  use  of 
(54),  we  arrive  at 

17  For  the  first  term,  Cij,  the  caveat  about  orbital  averaging  is  unimportant,  because  rj  bears  no  dependence 
upon  the  mean  anomaly. 
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dE 

r  .2  / 

,  15 

c  ■  -  ( 

,  121  2\1 

~d7  =  ~z 

1+2e) 

—  6  n  f]  e  1 

+  0(i2/Q)  +  0(enAt/Q).  (63) 


Thus  we  see  that  a  dissipative  tidal  torque  influences  the  spin,  while  causing  changes  also  in 
the  first  integral  E  (which  is  not  identical  to  the  actual  kinetic  energy).  Considering  equation 
(54)  and  noting  that  the  cosine  term  is  periodic,  we  see  that  weak  tidal  dissipation  will  cause 
a  change  in  i/2  over  time  scales  long  compared  to  P . 

For  nonsynchronous  rotation,  weak  dissipation  causes  both  a  slow  secular  change  and  a 
small  oscillation.  The  latter  arises  from  the  oscillating  part  of  t)  as  1 7  circulates  or  librates. 

Goldreich  (1966)  also  averages  over  the  period  P .  For  libration,  the  so-averaged  rate  of 
change  of  E  is 


( dE/dt  )p  =  -  2^~  U  +  y  e2J  +  0(e6)  +  0(i2/Q)  +  0(enAt/Q),  (64) 


where  the  positive  definite  quantity  W  is  introduced  via  (60).  The  negative  (dE/dt) p  for 
libration  means  that  the  maximum  i]2  at  cos  2rj  =  1  in  Eq.  (54)  is  decreasing  and  the  free 
libration  damps  with  time.  Libration  evolves  toward  synchronous  rotation  with  tj  =  0  and  rj 
equal  to  0  or  to  it . 

The  f]  term  causes  the  damping  of  the  libration  amplitude.  When  the  amplitude  is  suffi¬ 
ciently  small,  its  decrease  obeys  the  exponential  law  exp(—  DL  t)  with 


D 


L 


z 

2 C 


(65) 


As  we  emphasised  in  the  paragraph  after  equation  (61),  inside  the  1:1  resonance  the  time 
delay  At  showing  up  in  the  expression  for  Z  should  be  appropriate  to  the  libration  frequency 
Xlib  which  may  differ  greatly  from  the  usual  tidal  frequencies  for  nonsynchronous  spin.  This 
choice  of  At  will  influence  the  value  of  Z. 

For  circulation,  averaging  of  Eq.  (63)  over  one  orbital  period  leads  to 


( dE/dt )  p 


Z 

\ 

15  2\ 

2 

,  121  2\  ■  ' 

~P 

W  1 

1  +  T') 

—  12  jr  n  e2  1 

1  +  -jjr  e  J  sgn  (il)p 

+  0(e6)  +  0(i2/Q)  +  0(enAt/Q)  (66) 


where  we  have  used  definition  (59a)  for  P  and  definition  (60)  for  W. 

While  W  =  P  (rf)  p  is  positive  definite,  the  second  term  on  the  right-hand  side  of  Eq.  (66) 
can,  for  circulation,  have  either  sign.  In  the  case  when  its  sign  is  positive,  the  expression  (66) 
for  (dE/dt)  p  will  vanish  for  W  equal  to 

)  +  0(e6)  +  0(i2  /  Q)  +  0(enAt/Q).  (67) 

No  matter  whether  W  begins  its  evolution  with  an  initial  value  larger  or  smaller  than  WstaU , 
it  will  never  cross  W stall . 

Another  important  value  of  W  is  the  one  corresponding  to  the  boundary  between  libration 
and  circulation,  Wb-  Below  we  shall  obtain  its  value  and  shall  explain  that  for  Wb  <  Wstan 
there  exists  a  positive  value  of  (fj)  p  at  which  the  evolution  of  a  circulating  ( j]) p  should  stall. 


w stall  =  12  7T  n  e2  \ 


1  + 
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6.5  The  free-libration  frequency 


Exploring  librations,  we  start  out  with  the  small-amplitude  case.  Replacing  sin  2r]  with  2rj  in 
Eq.  (53),  we  write  down  the  frequency  for  small  librations  of  ry. 


Xlib- 


2jr 


B  -  A 


Ms, 


C  Ms, 
B  -  A 


+  Mp 
Ms, 


G 200 (e) 


1/2 


(68) 


In  many  realistic  situations,  the  condition  —  r  n  jj — 1  is  fulfilled,18  which 

ensures  that  Xub-max  ^  n-  Larger  amplitudes  increase  the  period  P  and  render  a  smaller 
frequency  Xm, ,  so  expression  (68)  gives  the  smallest  librating  period  and  the  largest  frequency 


Xlib-max  =  maxXlib- 


(69) 


The  linear  r)  term  in  (61)  will  alter  the  frequency,  but  for  slow  tide-caused  evolution  the 
correction  will  be  small,  so  (68)  still  will  serve  well  as  an  approximation  for  the  maximal 
frequency  of  libration.  On  all  these  grounds,  we  now  accept  that  for  both  small-amplitude 
and  large-amplitude  librations  the  inequality 


Xub  <  n  (7°) 

holds.  This  justifies,  a  posteriori,  the  assertion  made  after  (59b). 

In  the  literature  on  the  physical  libration  of  the  Moon,  the  expression  for  the  free-libration 
frequency  is  ubiquitous,  though  often  without  the  mass  factor.  Versions  of  the  expressions 
for  ribias  and  the  free-libration  damping  rate  appeared  in  Williams  et  al.  (2001).  That  paper, 
though,  did  not  rely  on  the  corrected  MacDonald  model  (constant  time  delay),  but  simply 
used  separate  Qs  for  libration  and  orbital  frequencies. 


6.6  The  boundary  between  circulation  and  libration 


The  boundary  between  circulation  and  the  resonance  zone  corresponds  to  a  value  W  =  W *. 
To  find  it,  we  combine  (54)  with  (56)  and  arrive  at 

1  7 

-  C  rf  =  Eh  (1  +  cos2tj),  (71a) 

which  is  the  same  as 

rf  =  A  cos2?;.  (71b) 

Insertion  of  the  resulting  expression  for  ?)  into  (60)  entails: 


Wb  =  8 


7r/2 


cos  rj  dr}  =  8  , 


=  4/7 


B  -  A 


Ms, 


Msec  4“  M  n 


G 200(e) 


1/2 


(72) 


comparison  whereof  to  (68)  yields  another  expression  for  the  boundary  value: 

Wb  =  4Xnb-max- 


18  This  condition  is  satisfied  safely  if  the  primary  is  a  planet  (like  Mercury)  or  a  large  satellite  (like  our  Moon). 
Its  fulfilment  is  not  guaranteed,  though,  for  small  satellites  (like  Phobos  or  Hyperion). 
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Fig.  1  Three  possible  scenarios 
of  evolution  of  {if)  p .  Arrows  on 
the  top  diagramme  illustrate  the 
evolution  of  circulating  (rj)  p 
toward  the  libration  boundary  for 
W rstaii  <  Wb-  Arrows  on  the  left 
of  the  bottom  diagramme  depict 
the  evolution  of  negative  [r])  p 
toward  the  libration  boundary. 
Arrows  in  the  midst  and  on  the 
right  of  the  bottom  diagramme 
show  the  evolution  of  positive 
(r])  p  toward  the  stall  point  for 
Wb  <  W,taU 


0 

Rate  <dr|/dt>p 


While  W  is  continuous  across  the  boundary,  P  has  a  logarithmic  singularity,  as  was  men¬ 
tioned  in  the  paragraph  after  equation  (56).  As  P  diverges,  the  evolution  rate  of  ( dE/dt  )p, 
given  by  (66) ,  vanishes  at  the  boundary.  Nonetheless  a  small  perturbation  allows  the  bound¬ 
ary  to  be  crossed — for  more  on  this,  see  the  two  paragraphs  after  formula  (22)  in  Goldreich 
(1966). 

6.7  Three  regimes 

If  Wstau  <  Wb,  then  there  is  no  stall  point  in  the  region  of  circulation  (Wb  <  W).  So  W 
evolves  towards  Wb ,  while  (fj)  p  of  either  sign  evolves  towards  zero  at  the  libration/circulation 
boundary.  Figure  la  illustrates  the  evolution  of  (fj)  P.  Goldreich  (1966)  comments  that  the 
boundary  will  be  crossed,  the  free  libration  will  damp,  and  the  rotation  will  evolve  toward 
the  synchronous  state.  The  synchronous  state  has  a  zero  r),  with  t]  biased  slightly  off19  of 
either  0  or  it .  The  figure  does  not  show  the  libration  region.20 

If  Wb  <  W^taii,  then  a  stall  point  exists  in  the  circulating  region  {Wb  <  W).  All  in  all, 
Fig.  lb  summarises  the  following  three  cases: 

A.  For  an  initially  negative  (i))p,  circulation  is  taking  place.  We  have  Wb  <  W,  but  the 
stall  point  is  never  located  on  the  negative  [f])  p  side  of  zero.  The  rate  (dE/dt )  p  is  neg¬ 
ative,  as  can  be  seen  from  (63).  We  then  have  a  slow  decrease  of  the  three  quantities:21 
E  — >  Eb,W  —>■  Wb,  and  (r/) p  —>■  0.  This  decrease  takes  the  system  to  the  circula- 
tion/libration  boundary.  Goldreich  (1966)  commented  that  it  is  ambiguous  whether  (fj)  p 
would  cross  zero  and  proceed  to  increase  or  whether  the  boundary  would  be  crossed 

19  The  value  jr  shows  up  because  of  the  factor  2  accompanying  rj  when  this  variable  enters  sin 2/)  in  the 
equation  for  the  torque  and  cos  2 rj  in  the  expression  E.  The  tiny  bias  off  rj  =  0  or  off  rj  =  tt  ,  given  by  (62), 
will  emerge  due  to  dissipation. 

20  A  figure  with  W  rather  than  (rj)  p  would  show  the  libration  region,  along  with  the  circulating  region.  On 
such  a  figure,  though,  positive  and  negative  values  of  (rj)  p  would  overlap,  so  that  the  one-sided  nature  of  the 
stall  would  not  be  evident. 

21  Be  mindful  that  E  remains  a  constant  over  a  cycle  of  P,  except  for  a  tiny  amount  of  dissipation  during  a 
cycle. 
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passing  into  the  libration  region  followed  by  damping  of  the  libration  and  by  evolution 
toward  the  synchronous  state. 

Be  mindful  that,  while  the  first  integral  E  is  decreasing  in  the  circulation  case,  the  actual 
kinetic  energy  of  rotation  is  increasing.  Indeed,  when  the  negative  [f]) p  is  evolving 
towards  zero,  the  spin  rate  6  is  growing,  and  so  is  the  rotational  energy  C62 /2. 

B.  For  an  initially  positive  (fj)  p  lying  between  the  circulation/libration  boundary  and  the 
stall  point,  i.e.,  for  Wb  <W<  WstaU,  the  rate  ( dE/ dt  )p  is  positive,  and  E  increases. 
The  quantities  W  ,  { f] 2)  p  and  { r])  p  will  evolve  to  larger  values  until  the  evolution  stalls 
as  W  approaches  WslaU  from  below. 

As  the  positive  r)  is  increasing,  so  are  the  spin  rate  9  and  the  rotational  energy.22 

C.  For  an  initially  positive  (fj) p  beyond  the  stall  point  we  have:  Wb  <  Wi(aH  <  W.Then 
{ dE/dt)p  is  negative  and  E  is  slowly  decreasing. 2l  Likewise,  W ,  (rj2) p  and  (ri) p 
evolve  to  lower  values  until  the  evolution  stalls  as  W  approaches  Wstall  from  above. 
The  decrease  of  the  positive  r]  renders  a  decrease  in  the  rotation  rate  9  and  thus  leads  to 
damping  of  the  kinetic  energy  of  rotation. 


Calculation  of  the  frequency  at  the  stall  point  ( r]stau  ) P  would  require  two  steps: 
first,  reversing  the  elliptic  integral  computation,  starting  with  Wstau ;  and  second,  deriving 
E/C  ,  Ps,aii  and  < i]stall  )P. 

At  the  stall  point,  formulae  (59a)  and  (60)  acquire  the  form  of 

Pstall  =  ,  l  71  .  and  W stall  =  Pstall  tftall)p>  (73> 

W stall  IP 

correspondingly.  In  combination  with  (67),  this  entails: 


Pstall )  P 


(74) 


whence  we  once  again  note  that  the  frequency  (f]stan  )  P  must  be  positive.  The  comparison 
of  WstaU  with  Wb  in  the  inequalities  mentioned  in  the  above  items  A,  B,  C  is  then  equivalent 
to  comparing  <ij2  q//  )P/{flstaii)p  with  i2/*)  Xlib-max  or  to  comparing24  3  n  e2  (1  +  e2/16) 

A  comment  on  Eq.  (74)  would  be  in  order.  Since  oscillations  of ;;  result  from  the  existence 
of  the  permanent  triaxiality,  evolution  of  becomes  smooth  in  the  A  =  B  limit.  Averaging 
becomes  unnecessary,  so  ^taII)p  becomes  simply  92nU,  while  (9stau) P  becomes  9stau- 
This  way,  in  the  oblate-body  case  considered  back  in  Sect.  5,  Eq.  (74)  acquires  the  form  of 
Pstall  =  6ne2(l  +  e2 / 16),  which  agrees  with  (44). 


6.8  Application  to  the  Moon 


When  Goldreich  (1966)  applied  his  inequality  expressions  to  the  Moon,  he  found  that  the 
stall  point  would  have  interrupted  evolution  of  rotation  from  faster  spin  to  synchronous 

22  In  the  first  integral  (54),  the  kinetic-energy-like  part  is  given  by  Cf/2 / 2  =  (C/2)(S  —  n)2  =  (C / 2)(92  — 
2n8  +  n2).  When  the  spin  accelerates,  the  additional  kinetic  energy  is  borrowed  from  the  orbital  motion. 
Because  of  the  dissipation,  the  overall  spin  +  orbit  energy  must  nevertheless  decrease.  (The  orbital  energy 
includes  a  kinetic  part  and  the  —GM/r  potential  term). 

22  Just  as  in  item  A,  here  E  stays  virtually  unchanged  over  P. 

24  The  latter  may  also  be  expressed  as  comparison  of  3 ne2  e2J  with  |^3  ^  ^  ^ 

M  Earth  1 1/2 

M Earth  “1“  M Moon  \ 
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rotation.  Here  we  have  repeated  his  study,  though  with  the  corrected  average  torque  (43) 
instead  of  (34).  For  —  ^  —  =  2.278  x  10-4  (Williams  and  Boggs  2012),  e  =  0.0549, 

and  -T-i — -  =  0.98785,  the  libration  period  turns  out  to  be  38  times  the  orbital 

M  Earth  +  MMoon 

period,  while  W stall  is  10  %  larger  than  Wb ■  This  renders  a  value  of  Wstaii  /  Wb  much  smaller 
than  the  one  found  by  Goldreich  (1966),  and  the  difference  is  mainly  due  to  our  use  of  the 
corrected  average  torque  (43).  Despite  the  so-different  value  of  Wstau/Wb ,  despinning  of 
the  Moon  would  still  be  interrupted  by  a  stall.  In  principle,  a  tidal  spin-up  scenario  remains 
an  option  too,  though  this  option  does  not  look  probable. 

Goldreich  (1966)  noted  that  the  lunar  orbital  eccentricity  is  changing.  To  make  the  stall 
point  disappear,  i.e.,  to  ensure  that  Wstaii  <  W*.  the  eccentricity  of  the  Moon  would  need 
to  be  less  than  95  %  of  its  present  value.  While  we  lack  data  on  the  ancient  evolution  of  the 
lunar  eccentricity,  the  modern  eccentricity  rate  of  about  2  x  10-11  year-1  has  been  reliably 
determined  through  analysis  of  the  Lunar  Laser  Ranging  data  (Williams  et  al.  2001 ;  Williams 
and  Boggs  2009).  For  the  measured  eccentricity  rate,  the  eccentricity  would  be  small  enough 
to  prevent  a  stall  prior  to  1 .4  x  1 08  year  ago.  The  Moon  has  clearly  been  a  satellite  of  the  Earth 
for  billions  of  years,  and  the  capture  into  the  synchronous  spin  state  should  have  occurred 
very  early  in  its  history. 

For  the  Moon  as  it  exists  today,  the  free  libration  in  longitude  has  a  2.9  year  period.  The 
damping  time  is  four  orders  of  magnitude  longer  (Williams  et  al.  2001),  so  evolution  of  this 
libration  is  slow.  Despite  the  damping  time  being  short  compared  to  the  lunar  age,  the  Moon 
has  a  small  free  libration  amplitude  of  1.3"  (Rambaux  and  Williams  2011).  There  has  been 
geologically  recent  stimulation,  probably  due  to  resonance  crossing  (Eckhardt  1993). 

The  ambiguity  of  evolution  of  the  negative  (17) p  past  the  circulation/libration  boundary 
deserves  a  comment.  For  Mercury,  Makarov  (20 1 2)  finds  that  a  more  realistic  tidal  dissipation 
model  than  the  corrected  MacDonald  torque  strongly  changes  computations  of  the  evolution 
of  planetary  spin  rate  near  the  3  :  2  and  higher  spin-orbit  resonances. 


7  Conclusions 

In  the  article  thus  far,  we  have  provided  a  detailed  explanation  of  how  the  empirical  MacDon¬ 
ald  model  can  be  derived  from  a  more  accurate  and  comprehensive  Darwin-Kaula  theory 
of  bodily  tides.  We  have  demonstrated  that  the  derivation  hinges  on  a  key  assertion  that  the 
quality  factor  Q  of  the  primary  should  be  inversely  proportional  to  the  tidal  frequency.  This 
crucial  circumstance  was  missed  by  MacDonald  (1964),  who  made  his  theory  inherently 
self-contradictory  by  setting  the  quality  factor  to  be  a  frequency-independent  constant. 

We  have  corrected  this  oversight  in  the  MacDonald  approach,  and  have  developed  an 
appropriate  correction  to  Goldreich’s  model  of  spin  dynamics  and  evolution  near  the  1:1 
spin-orbit  resonance.  Although  we  got  different  numbers,  qualitatively  the  main  conclusion 
by  Goldreich  (1966)  stays  unaltered:  when  an  oblate  body’s  spin  is  evolving  toward  the  reso¬ 
nance,  vanishing  of  the  average  tidal  torque  still  implies  rotation  slightly  faster  than  resonant, 
while  synchronicity  requires  a  small  compensating  torque.  For  a  triaxial  body,  the  picture 
gets  more  complex  due  to  the  emergence  of  a  triaxiality-caused  torque.  (While  the  oblate 
case  is  appropriate  for  spinning  gaseous  or  liquid  planets  and  moons,  the  triaxial  case  applies 
to  rocky  objects). 

Goldreich  (1966)  linked  the  possible  trapping  of  a  body  in  the  synchronous  state  during 
its  tidal  evolution  of  rotation  to  its  triaxiality.  In  light  of  the  correction  of  expression  (36)  to 
(44),  that  limiting  condition  between  the  triaxiality  and  e 2  must  change.  After  capture  into 
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the  synchronous  state,  the  tidal  torque  is  compensated  by  a  triaxiality  torque  by  aligning 
the  principal  axis  associated  with  the  smallest  moment  of  inertia  slightly  off  of  the  mean 
direction  to  the  external  body.  A  constant  is  thereby  introduced  into  the  physical  libration  in 
longitude. 

Setting  the  tidal  quality  factor  to  scale  as  inverse  frequency  is  incompatible  with  the 
actual  dissipative  properties  of  realistic  mantles  and  crusts.  Nevertheless,  after  the  afore- 
explained  correction  is  implemented,  both  the  corrected  MacDonald  description  of  tides  and 
the  dynamical  theory  based  thereon  remain  valuable  toy  models  capable  of  providing  a  good 
qualitative  handle  on  tidal  dynamics  over  not  too  long  timescales.  Specifically,  this  approach 
renders  a  simple  qualitative  description  of  the  interplay  between  the  tidal  torque  and  the 
triaxiality-caused  torque  exerted  on  a  body  near  the  1:1  spin-orbit  resonance. 

It  should  also  be  remembered  that  the  stability  of  pseudosynchronous  rotation  hinges  upon 
the  dissipation  model  employed.  Makarov  and  Efroimsky  (2012)  have  found  that  a  more  real¬ 
istic  tidal  dissipation  model  than  the  corrected  MacDonald  torque  makes  pseudosynchronous 
rotation  unstable.  Finally,  pseudosynchronism  becomes  impossible  when  the  triaxiality  is  too 
large  (see  subsection  6.7,  specifically  footnote  24). 
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Appendices 


Appendix  A:  The  tidal-torque  vector  and  its  components  in  spherical  coordinates 


A  secondary  with  spherical  coordinates  r*  =  (r*,  X*,  cf>*)  and  mass  M*ec  raises  a  tidal 
bulge  on  the  primary.  The  gravitational  attraction  between  the  tidal  bulge  and  a  secondary 
at  r  =  (r  ,  X  ,  0)  with  mass  Msec  causes  equal  but  opposite  torques  on  the  primary  and  the 
secondary.  For  the  external  tidal  potential  U  (r),  the  torque  components  depend  on  the  partial 
derivatives  of  the  potential  U  along  great  circle  arcs. 

To  calculate  these  expressions,  let  us  recall  some  basics.  The  torque  T  wherewith  the 
primary  is  acting  on  the  secondary  is  given  by  the  cross-product 

T  =  fxf,  (75) 


T  being  the  tidal  force  exerted  by  the  primary  on  the  secondary.  This  force  is  given  by 


T  =  -  AT  , 


dU(T) 

dr 


=  ~  M,. 


dU 

dr 


1 

+  - 


1  dU 


1  dU 


r  cos  0  dX 


~  77  .(76) 


r  30 


e^  ,  being  the  unit  vectors  of  a  spherical  coordinate  system  associated  with  the  pri¬ 


mary’s  equator  and  corotating  with  it.25  Insertion  of  (76)  into  (75)  results  in 


25  Were  we  using  the  polar  angle  ip  =  jr/2  —  <p  instead  of  the  latitude  <p,  the  right-handed  tri¬ 


ple  of  unit  vectors  would  be:  e 
+ 


ex  (“radial 

_  _ 1  3£7  a 

3  f  3r  T  r  dip  T  r  sin  ip  8/.  X ' 

As  we  are  employing  the  latitude,  the  right-handed  triple  changes  to  e 


<P  ' 

_  a  +  I 
r  r 


south  -  east”),  while  the  gradient  would  read  as 


e-  ,  (“radial  -  east  -  north”),  and 


the  gradient  becomes 


+ 


1  1 


8U 


r  cos  ip  dX  X 
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T  —  —  Msec 


( 


) 


(77) 


The  torque  r  wherewith  the  secondary  is  acting  on  the  primary  will  be  the  negative  of  (77): 


(78) 


Its  east  component,  the  one  aimed  along  ex,  is  parallel  to  the  primary’s  equatorial  plane. 
The  north  component,  aimed  along  ,  is  tangent  to  the  meridian.  In  our  paper  however  we 
employ  the  projection  of  the  torque  vector  onto  the  primary’s  spin  axis,  i.e.,  the  component 
r_  orthogonal  to  the  equator  plane.  This  component  is  cos  <p  times  the  north  component: 


dU 

al 


(79) 


This  torque  component,  Eq.  (6),  slows  down  the  rotation  rate  of  the  primary.  The  decelerating 
torque  acting  on  the  secondary  has  an  opposite  sign  and  is  rendered  by  (5). 

The  east  component  and  the  projection  of  the  north  component  onto  the  equator  plane 
act  to  change  the  orientation  of  the  primary’s  spin  axis,  a  subject  we  shall  not  pursue  in  this 
paper. 

Appendix  B:  Calculation  of  the  average  torque  within  the  corrected  MacDonald  model 

To  calculate  the  orbital  average  of  the  tidal  torque  acting  on  a  librating  secondary  obeying 
the  corrected  MacDonald  tidal  model,  substitute  (38)  and  (18)  into  (22),  and  then  choose 
a  =  —  1  (and  recall  that,  for  a  =  —  1,  the  integral  rheological  parameter  £  is  simply  the  time 
lag:  £  =  At).  An  equivalent  option  would  be  to  employ  (23)  directly.  This  will  lead  us  to 
the  following  expression  for  the  torque: 


r 


=  -\GM;eck 2\l£  \  9-v\  sgn(8-v)  +  0(i2/Q)+0(enAt/Q)  +  0(Q-3) 

v 


r 


r 


and  for  its  average  over  one  orbiting  cycle: 


3  GmIc  k2  £  d 


6 


,  3  GM;ec  k2£ 


R 


^j  +  0(i2/Q)+0(Q-3) 


+  O  (en  At  /  Q) 


(81a) 


+  e  cos  v)4  dv 
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|  3  GM]ec  h£ 


R 


(1  +e  cosv)6dv  +  0(r/Q)  +  0{Q^) 


0 


+  O  (en  At  /  Q) . 


(81b) 


Evaluation  of  the  integrals  is  straightforward  and  entails  (39-42). 
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